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Abstract 

It has been realized recently that in order to have a high level of 
maneuverability, supersonic delta wings should have a cross flow that is free 
of embedded shock waves. The conical cross flow sonic surface differs from 
that of plane transonic flow in many aspects. Well-known properties such as 
the monotone law are not true for conical cross flow sonic surfaces. Using a 
local analysis of the cross flow sonic line, relevant conditions for smooth 
cross flow are obtained. Using a technique to artificially construct a smooth 
sonic surface and an efficient numerical method to calculate the flow field, 
cones with smooth cross flow are obtained. 
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INTRODUCTION 


It is well known that the most suitable structure for a slender wing 
airplane has leading edge separation [1]. The resulting vortices are highly 
stable, and while contributing to the induced drag, also increase the lift 
because of the low pressure they induce. These vortices usually increase the 
lift-to-drag ratio for a delta wing. However, the modern supercruiser concept 
demands efficient supersonic cruise and high-level supersonic, as well as 
transonic, maneuverability [2]. In order to obtain this level of performance, 
the wings have to be relatively thick and should have transonic leading edges 
and attached flow [3]. A delta wing with transonic leading edge will, in 
general, have embedded cross flow shocks. To obtain attached flow, this 
configuration should have a cross flow that is free of embedded shocks or 
should contain only weak shock waves. Present study is tailored to find such 
configurations but is limited to conical wings and irrotational flow. Using a 
novel technique to construct smooth embedded sonic surfaces, shock-free cross 
flow is constructed for the first time for conical wings. We have also 
provided, based on a local analysis, the necessary geometric condition that 
should be satisfied by a configuration in order for a shock-free cross flow to 
exist. 


IRROTATIONAL CONICAL FLOWS 

Let {5°} be the surface coordinates of a unit sphere centered at the 

t v 

apex of the cone and g Q g be the corresponding metric tensor. Let v be 

the tangent velocity field (Figure 1) generated on the sphere by the 
mainstream flow. Irrotational assumption implies 
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where F is the conical potential. The total velocity q is given by 


q 2 - V° V a + F 2 . 


The continuity equation is [4] , 


(PV a ) |o + 2pF = 0, 


( 1 ) 


where II denotes surface covariant differentiation and 
and is given by the Bernoulli equation, 


P is the gas density 


P Y_1 = 1 +I|I M 2(l- q 2 ), 


( 2 ) 


Here M m is the free stream Mach number and y is the ratio of specific 
heats. Combining equations (1) and (2), we obtain the quasilinear form of the 
governing equation. 


r a3 V V f 2\ 

( g 2 — ^ V all 3 + ( 2 - M C ) F " °» (3) 

a 

where a is the speed of sound and M c is the cross flow Mach number 

( M c “ 1116 equation (3) changes from elliptic to hyperbolic type when 

M c increases through unity. The above behavior is similar to plane transonic 
flows and therefore can be utilized to develop an extremely efficient 
numerical method to calculate supersonic flows past arbitrary conical shapes 
at incidence [5] . 
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CROSS FLOW SONIC SURFACES 

Sonic bubbles that appear in conical cross flows differ from those of 

transonic flows in many aspects. In transonic flow the flow properties attain 
uniform state on the sonic line, whereas in conical flows the appearance of 
the radial velocity term in the Bernoulli equation causes the flow properties 
to vary along the sonic line. Plane transonic flow is governed by a 
homogeneous set of partial differential equations, and therefore, it is 
possible to obtain a linear problem through a Legendre transformation. This 
property also gives us the well-known Nikolskii-Taganov monotone rule [6] : if 

an observer moves along the sonic line, keeping the subsonic zone always to 
his left, then the stream vector will rotate in the clockwise direction. 
However, in conical flows, due to the inhomogeneous terms in the governing 

equations, it is not possible to obtain linear equations using a hodograph 
transformation, and also it is difficult to say anything definite about the 
streamline slope. Some properties of conical sonic surfaces have been worked 
out by Salas [7], and with the aid of some of his results, we will study some 
relevant aspects for shock-free flows. 

First we note an interesting behavior of the pressure at the point where 
a cross flow streamline exits the hyperbolic zone. Consider the Bernoulli 

equation 

-^-(q 2 + F 2 ) + — = constant; 

2 M c Y-l P 

taking the derivative in the s direction and substituting the adiabatic 

relation 

1 Y 
P - — 2 P ’ 

Y*T 

00 


we get 



9q 

_pq c^ + F) * 


Now, from the definition of the cross flow Mach number 


9M 
c 

9s 


= - [■ - M I s -] . 
a ds C d s 


Using this and the energy equations, we obtain the relation at the sonic 
surface for the streamwise pressure gradient, namely 



3M 

c 

3s 


+ f]. 


From this we see that since F is always positive, the pressure always 
decreases when the streamline enters the hyperbolic zone, which is similar to 
transonic flows. However, when it exits the hyerbolic zone, the pressure 
could either increase or decrease. This is in contrast to what happens in 
transonic flows, where the pressure always increases as the streamline exits 
the sonic surface because the radial velocity term F does not appear in the 
equation for 9p/9s. This possibility leads us to believe that a shock-free 
situation is more likely in conical flows. We recall here that in transonic 
flows the shock-free situation has been proved to be mathematically isolated. 
(It has been established experimentally [9], [10] that the neighboring flows 
have only weak shock waves.) A similar perturbation theory for shock-free 
conical flows does not yet exist and shock-free solutions have not previously 


been shown to exist. 
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Next, we will observe a unique feature of the cross flow sonic line. We 
will first obtain a description for the sonic line by expanding the cross flow 
velocity and the speed of sound and equating them. Thus we write 


* 

a = a + 



q 


c 



• • • 


and also, from the energy equation, 


(ft* - - (V'jfe*- 


2 ^3n 


Thus, the equation of the line q = a is 


(— 1 . . = - r_L. 

'■ds' sonic line V Y + 1 


3M * 

{—) 

* l 3s j 


3q * 
(— 1 ) 
l 3n J 


This relation indicates that the angle at which the cross flow sonic line 

3q * 

meets the body is determined by the sign of (^— ) • 

This means that if this quantity is negative, the sonic line could meet 
the body at an obtuse angle. Salas incorrectly predicts this possibility for 
a circular cone and then uses an argument based on pressure to speculate on 
the possibility of a shock wave. However, one could rule out the existence of 
smooth flow with this kind of a bubble, using the following argument. 
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At the sonic line, the characteristics are normal to the streamlines, and 
therefore, from Figure 2, it is clear that for a sonic line of this kind, 
characteristics of the same family will intersect. Thus, if through the above 
analysis, we arrive at a sonic line of this type, then for this cone a shock- 
free cross flow is impossible. In fact, one may work out rules to identify 
conical shapes for which shock-free cross flow is impossible. For 
convenience, we will work in the stereographically projected plane. If one 
chooses a Cartesian system (x,y) in this plane, then the metric tensor will 
be 


Sag 


where 


J = 


, . 2 . 2 
1 + x + y 


If a denotes the iclination of the streamline with 


in = 

dx 


tan a. 


then the irrotationality condition becomes 


3q c r_9a _1 9 j -| 

9n “ V 9s “ J 3n 


We have already noted that when 



9n 


< 


0 


at the body, shock-free cross flow on a general cone is impossible, and hence 
cones with 


1 3 J 

R " Tn 


< 0 


shock-free cross flow is impossible. 
We note that 


3 1 A 3 ^ 1 . * 3 

= -7 cos a — + rr sin a tt 

3s J 3 x J 3 


3 1 . A 3 A 1 A 3 

*— = — r sin a tt— + — cos a 

3n J 3x J 3 


, 2 _ 2 ( 2 A , 2v 

ds = J (dx + dy ) , 


and R is the radius of curvature of the image of the cone in the 

stereographical ly projected plane. Thus, cones that satisfy the above rule 
will not have a shock-free cross flow. 

We will now show that this condition is not satisfied on a circular 
cone. Consider a circular cone of half angle • Then the radius of the 

circle on the stereographi cally projected plane is 

„ *c 

a = tan — , 


since 


i|> c < 90°, a < 1. 


Thus on the body 


1 + a 


J 
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and 


_ q c f 1 . ^ 

8n J U Ja ^ 


> 0 . 


A DESIGN METHOD TO OBTAIN CONICAL WINGS WITH SHOCK-FREE CROSS FLOW 

It is possible to devise a direct approach to search for shock-free 
configurations. Suppose we consider a test case with shocked cross flow and 
change the expression for the density (equations) inside the cross flow bubble 
in such a way that the resulting partial differential equation is elliptic; 
then the ellipt ic-to-elliptic transition will result in a smooth sonic 
surface. 

We may accept the solution outside this surface (which includes the bow 
shock wave and part of the cone) and use the flow properties on this surface 
to solve the Cauchy problem for the actual gas law (for the actual gas law, 
the governing equation is hyperbolic) to obtain the new body shape inside this 
surface. However, one should note that in this method there is no guarantee 
that a certain gas law will provide Cauchy data that will provide a smooth 
flow up to the body. This method is still preferable because it is direct and 
only part of the configuration is being modified. This method has been 
successfully introduced to transonic flows by Sobieczky [11]. The application 
of this fictitious gas method is not straightforward for conical flows for the 
following reasons. 

In plane transonic flow, a sonic bubble is defined by the statement 
q > a , where a is the speed of sound at the sonic condition and is a 
known constant. In conical flows, however, a varies and thus is an unknown. 
This difficulty can be eliminated by first computing the actual speed of sound 
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and, whenever it is less than the cross flow velocity, replacing it by the 
fictitious speed of sound. When we change the gas law, we must take care to 
preserve mass conservation at the sonic line. In conical flows, because the 
density varies along the sonic line, the gas law should be chosen to give 
continuous density across the sonic line. Let us look at the simplest way to 
meet these requirements. The energy equation is 





and at the sonic conditions 


M 


2 _ y- 1 
a * 


1 + 


Y-1 


M 2 (1-F 2 ) 

oo "k 


P£) 


Thus, if we use the gas law of the form 



1 


+ 


Y-1 

2 


m 2 (i-f 2 ) 


( 





then the flow properties would be continuous across the sonic line and the 
resulting partial differential equation will have the form 


g 


cxg 3 F 


3 ^ 


= 0, 


and be elliptic. This could be called the incompressible analogue for conical 
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flows. In changing the gas law inside the sonic surface, one needs to make 
sure that the gas law near the bow shock wave is correct. 

This was done by first solving the real problem to convergence, so that 
the bow shock and the cross flow sonic bubble were well developed, and then 
using this as the initial condition to solve the problem with the fictitious 
gas law inside the sonic surface. Inside the bubble the artificial viscosity 
[5] should be switched off and the iteration scheme should be specialized to 
the case of a + We will now consider the Cauchy problem for the cross 

flow bubble. The domain Q is first mapped to a rectangle, using an '.'onion 
peel transformation as shown in Figure 3. Since this transformation is 
singular at the end points, it is decided to use the covariant velocities in 
the spherical polar coordinate system (0 ,xp) as dependent variables. 

Let 

U = F q and V = F, : 

0 y 7 

then the governing equation is 

AU 0 + B(U. + V.) + C V, + D = 0 
0 0 

and the irrotat ionality gives 


where 



sin a sin i}> 
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B 


2 

a 


UV 


sin 
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C = 1 



2 

a 


and 

2 

D - (2 - M 2 )F + V cot i|>( 1 + - 5 — 2-= — ) . 

c a sin ^ 


The derivatives of covariant velocities are now transformed from (0, 
(C,n) plane. Thus we get 


U n 


P + 


Here 




The Cauchy problem is then solved by a predictor-corrector method* We 


f* = f° + Anfve^ + vi|> n )° 


u 


0 , 0 — -K) -*0\ 

U + An(p ( 6 ^ U ) + Q J 


F + = y[F* + F° + An(U0 n + V* )*] 


to 


write 


IT = + U + Ati ( p u ) + Q )] 
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ON THE WELL-POSEDNESS OF THE MARCHING PROBLEM 

The marching problem would be well posed in the absence of any kind of 
singular points or limit lines. Conical equations do allow the cross flow 
streamline to meet the sonic line at right angle. However, at this point, the 
characteristics are tangent to the sonic line and the marching is not well 
posed. This situation occurs for two kinds of sonic lines, as shown in Figure 
4. Only type (B) is relevant for our design. From the earlier arguments we 
see that for this kind of bubble, the same family of characteristics will 
intersect, and therefore if the fictitious gas chosen leads to a bubble of 
this shape, it should be discarded. 

A sonic line that does not have such a singular point on it could still 
lead to a limit line (as shown in Figure 5) in the subsequent marching. If 
this limit line occurs above the body, then this sonic line should again be 
discarded. 

The Jacobian's 


J„ = 9 \J> - 9 tj>_ 

g C n n 


and 



- 2B "e 



5 


respectively, signal grid singularities and limit lines by changing their 
signs. 
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RESULTS 

A simple test case is given to validate the method. A circular cone of 
10° half angle at an angle of attack of 20° at Mach 2 is considered. Figures 
6 and 7, respectively, show the sonic surfaces and pressure distributions for 
the original as well as the modified cone. Figure 8 shows the required 
surface modification of the cone. In the actual design process the fictitious 
gas method should be combined with a change in camber to produce the required 
performance. 
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Figure 1 


Figure 2 


Figure 5 


Figure 6 


Figure 7 


Figure Captions 

The classical tangent field representation of conical flows on unit 
sphere. 

Types of sonic lines that terminate in the body. 

Limit line. 

Bow shock position and the cross flow sonic line for a 10° circular 
cone at 20° angle of attack at M =2. 

original cone, • modified cone 

Surface pressure distribution on the circular cone. 

□ original cone, o modified cone 


Figure 8. 


Surface modification of the 10° cone. 
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Figure 5 
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